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We measure the statistical distribution of the photonic local density of states in the near field of a 
semi-continuous gold film. By varying the distance between the measurement plane and the film, we show 
that near-field confined modes play a major role in the width of the distribution. Numerical simulations in 
good agreement with experiments allow us to point out the influence of non-radiative decay channels at short 
distance. © 2012 Optical Society of America 
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The local density of optical states (LDOS) is a quan- 
tity of fundamental interest for the understanding and 
control of light-matter interaction in structured envi- 
ronments [1-3]. Plasmonic structures are expected to 
produce substantial changes in the optical modes both 
in space and frequency, offering degrees of freedom for 
LDOS engineering on the nanometer scale at visible and 
near infrared wavelengths. The LDOS can be probed us- 
ing fluorescent emitters [4-6] or thermal radiation [7, 8] . 
Qualitative LDOS mapping can also be performed in 
scattering experiments using angularly-integrated inco- 
herent illumination or detection [9,10]. 

The LDOS p(ro,w) at position ro and frequency a; is 
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where G is the dyadic Green function, such that the elec- 
tric field radiated at point r by a point electric dipole 
p located at point rg reads E(r) = /iow^G(r, Tq, a;)p. 
This LDOS is relevant for the spontaneous emission by 
an electric dipole emitter with position rg and emission 
frequency w, and without defined orientation [3]. The de- 
cay rate averaged over the orientations of the transition 
dipole p reads: 
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Measuring the fluorescence lifetime t = l/F is a di- 
rect way to probe the LDOS. In plasmonics, the seminal 
work of Drexhage demonstrated substantial changes in 
the lifetime when approaching a flat silver surface [1,11]- 
More generally, the distance dependence of fluorescence 
signals in the near held of flat surfaces [12], isolated 
nanoparticles [13, 14] or nanoparticlcs arrays [15] reveals 
the interplay between radiative and non-radiative pro- 
cesses. Understanding this interplay is of high interest 
for the design of efficient light sources (radiative pro- 
cesses), or for the control of fluorescence quenching and 
nanoscale energy transfer (non-radiative processes). 

In this Letter, we study the distance dependence 
of the LDOS in the near field of a disordered semi- 



continuous gold film. Such fllms are known to sup- 
port unconventional plasmon modes, that, in particu- 
lar, can localize optical energy in subwavclcngth areas 
(hot spots) [5,16-19]. We measure statistical distribu- 
tions of the LDOS at different distances in the near field, 
and compare the experimental data to numerical simu- 
lations. This allows us to assess the role of radiative and 
non-radiative channels. 

The experimental setup and the sample geometry are 
sketched in the inset in Fig. 1. The sample prepara- 
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Fig. 1. Distributions of the decay rate (normalized by 
its averaged value on a glass cover-slide) for differ- 
ent distances d between the center of the fluorescent 
nanosources and the film. Inset: sketch of the setup. 

tion consists in the following steps. First, gold is evap- 
orated on a microscope cover-slide, in order to obtain 
a semi-continuous film around the percolation thresh- 
old. Second, a Si02 spacer is evaporated on the gold 
film. Third, fluorescent nanosources (polystyrene beads, 
Invitrogen Fluospheres Red 580/605, diameter 25 nm) 
are spin-coated on the top of the silica spacer. Sev- 
eral samples with different spacer thicknesses ranging 
from 20 to 80 nm have been produced. The experimen- 
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tal setup used to measure the fluorescence lifetime of 
the nanosources has been described previously [5]. We 
briefly recall the main features. The nanosources are 
individually illuminated through an oil immersion ob- 
jective mounted on an inverted microscope. The excita- 
tion is performed at a repetition rate of 10 MHz by a 
picosecond super-continuum laser (Fianium SC450) fil- 
tered at the excitation wavelength 560 nm, with an in- 
tensity ~ 10^ W/cm^. Fluorescence photons are col- 
lected through the same objective. They are separated 
from excitation photons by a dichroic mirror followed by 
a set of filters centered around A = 607 nm, with a width 
of 70 nm, and are detected by an avalanche pliotodiode 
(Perkin Elmer SPCM-AQR-15). Lifetime measurements 
are performed by time-correlated single photon counting 
(PicoQuant TimeHarp200). 

Figure 1 displays the distributions of the decay rate at 
several distances d between the center of the nanosources 
and the gold film, obtained from lifetime measurements 
on 30 nanosources. One clearly sees that both the mean 
value and the width of the distribution are affected by 
the proximity of the disordered gold film. The behavior 
of the averaged value can be qualitatively understood 
replacing the semi-continuous film by an effective ho- 
mogeneous film. The averaged LDOS close to this ab- 
sorbing homogeneous film is expected to be larger than 
that on the glass cover-slide (our reference) due to non- 
radiative channels, and to increase at short distance [1]. 
The broadening of the decay rate distribution when the 
distance to the film decreases is more interesting. The 
disordered nature of the film induces LDOS fluctuations, 
that are known to become substantial in the presence of 
near-field interactions able to generate large field varia- 
tions on scales smaller than the wavelength [5,20]. An 
important feature in our data is the substantial change in 
the shape and width of the statistical distribution with 
the distance to the film, in the range c? ~ 30 — 90 nm. 
This is due to spatial filtering of optical modes later- 
ally confined on scales below the wavelength. Indeed, the 
field distribution in a plane at a distance d is exponen- 
tially filtered in Fourier space by a factor ^ exp{—Kd) 
compared to the distribution at d = nm, with K 
the spatial frequency in the transverse direction (par- 
allel to the film plane). This is a feature of near- field 
optics [3] . From this simple observation, an order of mag- 
nitude of the lateral confinement ^ of the field can be ex- 
tracted. Since field variations giving rise to substantial 
fluctuations of the decay rate strongly attenuate between 
d — 32.5 nm and d = 52.5 nm, the attenuation length 
can be estimated to be 1/K ~ 10 nm. One can deduce 
^ ~ 2'k/K ~ 60 nm as a typical size of hot spots, in 
agreement with orders of magnitude found by near-field 
optical microscopy [18,19]. 

In order to have a better insight into the physical phe- 
nomena driving the decay rate distributions, we have 
performed numerical simulations. We have generated 
semi-continuous gold films using a Kinetic Monte-Carlo 
algorithm [21] that mimics the evaporation and depo- 



sition process. The principle is to let 5 nm gold parti- 
cles, interacting via a rescaled atomic potential taken 
from [22], diffuse on a 2D square matrix until a stable 
configuration is reached. A typical realization of a film 
with feature comparable to that used experimentally is 
shown in the inset in Fig. 2. Numerically generated films 
exhibit fractal clusters close to the percolation threshold, 
as observed in TEM images obtained on real films [5]. 
To calculate the electric field at any position r above the 
film, under illumination by a point electric dipole p lo- 
cated at position Tq, we use a volume integral method 
based on the solution of the Lippmann-Schwinger equa- 
tion: 

E(r) = Eo(r) + ^^[e(c^)-l] / Go(r, r', c.)E(r')dV 
c Jv 

(3) 

where Gg is the vacuum dyadic Green function, e(w) the 
dielectric function of gold, V the volume occupied by 
the gold film and Eo(r) = /.Jow^Go(r, ro, w)p is the inci- 
dent field. The numerical computation is done by a mo- 
ment method, discretizing V into 2.5 nm cells. The Green 
function Gg is integrated over the cell to improve conver- 
gence [23]. Once the electric field is known, the full Green 
function G is deduced from E(r) = /.Joa;^G(r, Tq, u;)p, 
and the decay rate is computed from Eqs. (1) and (2). 

Calculated distributions of decay rates are shown in 
Fig. 2, for the same distances d as in Fig. 1. The 
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Fig. 2. Calculated distributions of the normalized decay 
rate F/Fo (Fq is the value in vacuum) for different dis- 
tances d, as in Fig. 1. A = 605 nm. Inset: Film used in 
the simulation (black color corresponds to gold). 

nanosources are modeled by orientationally averaged 
electric dipoles distributed inside spheres with diameter 
25 nm. The metal film and the nanosources are immersed 
in a homogeneous glass background that models the glass 
substrate and the silica spacer. Experimental and numer- 
ical data show a very good qualitative agreement for the 
trends of both the mean and the width of the distribu- 
tions. An interesting aspect of numerical simulations is 
the possibility of computing separately the radiative and 
non-radiative decay rates Fr and Fnr. The radiative de- 
cay rate is proportional to the power radiated in the far 
field, while the non-radiative decay rate is proportional 
to the power absorbed in the sample [13,24]. The mean 
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value and variance of the distributions of T/Tq, Tb,/To 
and rNR/ro are represented in Fig. 3 versus the distance 
d. The increase of (T) at short distance is mainly due 
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Fig. 3. Mean value (left) and variance (right) of the dis- 
tribution of F/Fq (black sohd line), Fnr/Fq (blue dotted 
line) and Fr/Fq (red dashed line) versus the distance d. 

to the increase of (Fnr). This can be understood recall- 
ing the effective homogeneous film approach. Indeed, one 
does not expect significant variations of (Fr) on such a 
distance range, since variations in the radiative rate are 
due to interference effects that build up on a length scale 
of the order of A/2 ^ d. Interestingly, the broadening 
at short distance is also driven by non-radiative decay 
channels. The confined near-field variations responsible 
for large near-field fluctuations of the LDOS are chiefly 
generated by non-radiative modes. 

The numerical calculation also allows us to compute 
spatial maps of both Fr and Fnr on a specific sample, as 
shown in Fig. 4. The similarity between the spatial dis- 
tribution of F and Fnr is striking, confirming the influ- 
ence of non-radiative decay channels on the LDOS vari- 
ations. The spatial distribution of Fr and Fnr do not 
coincide, showing that some specific points on the sample 
might give rise preferentially to a coupling to radiative 
or non-radiative channels. Finally, we point out that the 
spatial extent ^ of bright LDOS spots is in agreement 
with the qualitative estimate based on the experimental 
data. 




Fig. 4. Maps of the normalized decay rates for a single 
orientationally averaged dipole at a distance d = 32.5 nm 
(a. b. c) and d = 52.5 nm (d, e, f) above the film shown 
in Fig. 2. (a, d): F/Fq, (b, e): Fnr/Fq, (c, f): Fr/Fq. 

In summary, we have studied the distance depen- 
dence of the statistical distributions of the decay rate 
(LDOS) in the near field of a disordered semi-continuous 
gold film. Qualitative analyses and numerical simulations 



have allowed us to put forward the influence of confined 
near-field modes, and to assess the role of non-radiative 
channels at short distance. This work illustrates the po- 
tential of LDOS statistics in the study of optical modes 
in nanostructured disordered materials. 
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